Effects of Transverse Magnetic Anisotropy on Current-Induced Spin Switching 
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Spin-polarized transport through bistable magnetic adatoms or single-molecule magnets (SMMs), 
which exhibit both uniaxial and transverse magnetic anisotropy, is considered theoretically. The main 
focus is on the impact of transverse anisotropy on transport characteristics and the adatom's/SMM's 
spin. In particular, we analyze the role of quantum tunneling of magnetization (QTM) in the mecha- 
nism of the current-induced spin switching, and show that the QTM phenomenon becomes revealed 
as resonant peaks in the average values of the molecule's spin and in the charge current. These 
features appear at some resonant fields and are observable when at least one of the electrodes is 
ferromagnetic. We also show that the conductance generally depends on the relative orientation of 
the average adatom's/SMM's spin and electrode's magnetic moment. This spin-valve like magne- 
toresistance effect can be used to control spin switching of the adatom's/SMM's spin. 



Experiments on electronic transport through individ- 
ual atoms/molecules are at the forefront of the search for 
novel nanoelectronics and information processing tech- 
nologies. In this context, very prospective are magnetic 
atoms p~H3] and single-molecule magnets (SMMs) (4j-[8] 
with a large spin S > \. If properly deposited onto a 
substrate, these quantum systems can acquire (in the 
case of atoms) [9] or retain (in the case of SMMs) fTU] 
their intrinsic magnetic anisotropy - a property respon- 
sible for magnetic bistability. Especially attracting is 
the idea of incorporating magnetic adatoms/SMMs into 
spintronic devices [TT], with the objective to use spin- 
polarized currents for manipulation of their magnetic mo- 
ments |12|[T3] . Actually, the feasibility of this concept has 
already been experimentally proven for Mn adatoms [3]. 
One of the key conditions for successful applications is a 
sufficiently large uniaxial magnetic anisotropy constant 
D. Therefore, some efforts have been undertaken in order 
to synthesize new molecules with large D or to find other 
ways of anisotropy enhancement. It has been also demon- 
strated that magnetic anisotropy of an adatom/SMM can 
be systematically tuned, albeit in a limited range, e.g., 
by the environment adjustment |14j . an external electric 
field |6J , or mechanical stretching of a molecule [7] . 

Apart from the uniaxial magnetic anisotropy un- 
derlying the magnetic bistability, adatoms and SMMs 
usually possess also the transverse component of the 
anisotropy [4J. If the latter component is sufficiently 
large, it may lead to additional quantum effects, like oscil- 
lations due to the geometric Berry phase |15| or quantum 
tunneling of magnetization (QTM) |16[[T7]. Although the 
role of QTM in electronic transport has been studied 
extensively for normal electrodes [18H20] , much less is 
known how it affects the spin-polarized transport |21| . 
Since QTM allows for the underbarrier transitions be- 
tween the states on the opposite sides of the energy 
barrier, it may serve as an additional dephasing mech- 
anism, and thus impede the control of spin state by spin- 



polarized currents. 

In this Letter we address the mechanism of current- 
induced spin switching in the presence of transverse 
anisotropy. We show that the conductance reveals peaks 
at voltages where the thermal transition rates directly 
between degenerate states of lowest energy are equal to 
the rate of transitions to the first excited state. More- 
over, the transverse anisotropy significantly modifies the 
current-induced spin switching at some resonant fields, 
where the QTM phenomenon leads to resonant peaks in 
the field dependence of the average value of spin and 
charge current. These resonances are well pronounced in 
the field dependence of the derivative of current with re- 
spect to the magnetic field, and appear only when at 
least one electrode is ferromagnetic. We also show that 
the conductance generally depends on the relative ori- 
entation of the magnetic moments of the electrode and 
adatom/SMM. This effect follows from the interference of 
direct and indirect spin-conserving tunneling processes. 
In addition, a significant bias reversal asymmetry appears 
then in the transport characteristics. 

Model. — Key features of magnetic adatoms and SMMs 
are captured by the giant-spin Hamiltonian [1], 

H s = -DS 2 z + j(Sl + S 2 _) + S-B, (1) 

where the first and second terms stand for the uni- 
axial and transverse magnetic anisotropy, respectively, 
while the last term represents the Zeeman interaction, 
with B = (B x , By, B z ) denoting an external magnetic 
field measured in energy units. Since we are interested 
here in systems with an energy barrier for spin switch- 
ing, we assume D > 0. Without losing generality, we 
also assume positive perpendicular anisotropy constant, 
E > 0, and s? E/D < § 0]. When E ^ 0, each 
of the 25 + 1 eigenstates \x) of the Hamiltonian (JTJ, 
T~is\x) = Ex\x)i i s a linear combination of the eigenstates 
| to) of the S z component. We label the states \\) with a 
subscript to, |%) — > \Xm)i which [as well as the numbers 
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modeled by the Appelbaum Hamiltonian [121 H51 1221 125] . 
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Figure 1. (Color online) (a) Schematic depiction of the 
system under consideration, (b) Examples of inelastic elec- 
tron tunneling processes due to scattering of electron spin 
on the adatom's/SMM's spin, (c) Energy spectrum of the 
adatom/SMM for S — 5/2 in the absence (left) and pres- 
ence (right) of an external magnetic field along the easy axis, 
B = (0,0, B z ). Bottom panel shows the differential conduc- 
tance as a function of bias voltage in the case of nonmagnetic 
electrodes (P t — P s — 0) for: (d) selected values of the uniax- 
ial magnetic anisotropy constants D and absence of transverse 
anisotropy, E = 0; (e) several values of E for a given D. Re- 
maining parameters: T<j = 0.1 eV, p* — p s = 0.5 eV _1 [thus 
Go « 0.025 a = 1 and 2r?< = tj, = 1. 



in Fig. [TJc)] corresponds to the S z component of highest 
weight in the state \xm), i- e - \Xm) = \ m ) f° r T— >0. When 
B = (0, 0, B z ), the eigenstates \Xm) can be written as 
\Xm) = Efc=o,± 2 ( TO + k \Xm) \m + k) for m = -S, . . . , S, 
where (m\xm) is the amplitude of the state |m) in the 
eigenstate \xm)- Thus, the transverse anisotropy leads 
to mixing of the states |m), and therefore enables the 
QTM [17J. In the following, the index m shall be used 
only when necessary to avoid any confusion. 

We consider an experimental configuration including 
tip of the scanning tunneling microscope (STM), and a 
substrate [which plays the role of second electrode] at 
which the adatom/SMM is deposited, see Fig.^a). Both, 
STM tip and substrate are characterized by noninteract- 
ing itinerant electrons, Hei = X) gk(r £ L a ko- a L ((? = * for 
the STM tip, and q = s for the substrate), with the energy 
dispersion e£ , and a£ (a^cr) being the relevant creation 
(annihilation) operators (k is a wave vector, and a is the 
electron spin index). In general, both electrodes can be 
magnetic, with spin-dependent density of states (DOS) 
p q a at the Fermi level. By introducing the spin polariza- 
tion coefficient, P q = (p^ — pj)/(p^ + p-p, the DOS can be 

parameterized as: p!Lj> = ^-(1 ± P q ) with p q = p^ + pj. 
Electron tunneling processes in the STM geometry are 



Ht = 



EH 



ka a k' 



q'0 



qkk'a 



(2) 



with q to be understood as s = t and i = s. Electrons 
can tunnel either directly between the two electrodes [the 
first term Eq. (pi], or during the tunneling event they can 
interact with the adatom/SMM via exchange coupling 
[the second term in Eq. pj]. The former processes are de- 
scribed by the tunneling parameter Td, whereas the latter 
ones by the exchange parameter J qq i , see Fig.[l];b). Both 
Td and J qq i are assumed to be real, isotropic, and inde- 
pendent of energy and the electrodes' spin-polarization. 
In the following discussion, we write J qq i — Jn q n q i, with 
ij q denoting the dimensionless scaling factor of the cou- 
pling between the adatom/SMM and the gth electrode 
(we fix rj s = 1). We also relate the parameters J and Td 
as J = aTd. Thus, Td will serve as the key, experimentally 
relevant parameter [23 , and a describes relation between 
the direct tunneling processes and those involving spin 
scattering of conduction electrons. For simplicity, elec- 
trodes' magnetic moments are assumed to be collinear 
with the easy axis of the adatom/SMM. 

Method. — In the weak coupling regime, transport char- 
acteristics can be derived using the approach based on 
a master equation. The charge current (e > 0) flowing 
between the STM tip and the substrate is then given 

by I = e E XX ' v x{lxx' ~ 7 xx'}' where p x is the proba- 
bility of finding the SMM/adatom in the spin state |x) 

(= \xm)), and 7**, stands for the transition rate between 
the states |x) and \x') (= \Xm')) associated with electron 
tunneling between the electrodes q and q' . 

For the sake of analytical clarity, we decompose the 
total current into two parts, I = / e l + -An, where J e i rep- 
resents the contribution due to elastic electron tunnel- 
ing processes [the spin state \x) remains unchanged] and 
/i n is the inelastic term [with transitions between dif- 
ferent states \x) and |x')]. In the second order approxi- 
mation with respect to the electrode-SMM/adatom cou- 
pling, these two current components take the form 

^ = {1 + P t P s + 2a(S z )(P t + P S )}V + J2^xx> (3) 
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In the equations above, Go 
and §*, x = (x'|5 K |%) for k = Z,±. Accordingly, (S z ) = 

Y.x^xKx- In addition- A eV = E x ~ E x' + eV ; where 
eV = fit — Ms stands for the difference in electrochem- 
ical potentials of the tip (fj, t ) and substrate (fj. s ), while 
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) with ((x) = x/{l 



exp(— ^/{/cbT 1 })} and T denoting temperature. 

In order to evaluate the current from Eqs. ^ and 
Q, we need the probabilities P x . These can be ob- 
tained from the set of stationary master equations, V x : 

{^x'7x'x ~ ^x7xx'| = 0; where the golden rule tran- 
sition rates j xx > = V 7w for X^x' are g iven b Y 
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Nonmagnetic electrodes. — Consider a model system of 
spin 5=5/2, connected to nonmagnetic tip and substrate, 
and characterized by typical parameters observed in ex- 
periments, see the caption of Fig. [T] For vanishingly small 
transverse magnetic anisotropy (E = 0) and |B| = 0, the 
Hamiltonian (JTJ is diagonal in the basis of the eigen- 
states of S z . As a result, Eqs. (pH) and Q simplify signif- 
icantly [H]: /| =0 /Go = {1 + &{S 2 Z )}V and lF=°/G = 

f e T. m V m Y.^ ± AM^)\ 2 Z-(Kv m+X ) with A\(m) = 
[S{S+ 1) -m(m + A)] 1 / 2 . For L> > and low T, the 
system occupies then with equal probabilities each of the 
metastable ground states |± S). When bias voltage \V\ 
increases, initially only the elastic tunneling processes 
contribute to transport, i.e. I^ =0 ^ and /^ =0 ~ 0. 
When \V\ becomes of the order of the threshold value 
V thl = D(2S - 1) = AD (for S = 5/2), see Fig. [jjc), 
the inelastic processes become activated and IP =0 ^ 0, 
which appears as a characteristic step in the differential 
conductance, Fig. [ljd,e). This feature is typical of the 
systems exhibiting easy-axis magnetic anisotropy |5J [5] . 

For E > and a half-integer spin S, the eigenstates 
\Xm) are twofold degenerate (Kramers' doublets) and 
form two uncoupled sets {|x±S=F2fc)} fc=0>1 s-i/2 PS, 
schematically distinguished by different colors in 
Fig.[IJc). Mixing of the states corresponding to different 
values of m is also revealed in the expressions for current, 

\Y, m \( m \Xk)\ 2 m\ 2 }V and 
{\J2m(xi\m)(m\xk)m\ + 

2-i 
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The transverse anisotropy manifests as several new fea- 
tures in transport characteristics. First, for \V\ < I4hr it 
leads to reduction of I e \ [cf. the second terms of /Jf =0 and 
I^i given above] , which appears as a decreased value of 
dl/dV, Fig. [TJe) . Second, the transverse anisotropy ef- 
fectively gives rise to a small reduction of T4hr- Actually, 
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Figure 2. (Color online) Dependence of differential conduc- 
tance (a) and the average value of the spin's 2th component 
(b) on bias voltage and transverse magnetic anisotropy for 
D = 50 /^eV and Pt = 0.5. Solid lines in (c) and (d) rep- 
resent cross-sections of (a) and (b), respectively, for selected 
values of E, while the dashed line corresponds to E/D = 
and Pt = 0. Other parameters as in Fig. [I] 



at equilibrium, V ~ 0, the spin can directly oscillate be- 
tween the two ground states |x±s) as 7x-sxs = 7xsx-s (x 
kftT |24) . which is in contrast to 7_ss = 7S-S = for 
E = 0. Such 'underbarrier' oscillations will dominate un- 
til they are surpassed by transitions to the first excited 
spin dublet, which occurs at \V\ ~ V t hi(E). It's worth 
emphasizing that the effect stems entirely from thermal 
fluctuations, so that the Kramers' degeneracy is not af- 
fected. Finally, for significantly large transverse magnetic 
anisotropy the second step in the differential conductance 
appears at V^. « 4D(S-1) = 6D (for S = 5/2), when di- 
rect spin excitations to the second excited Kramers' dou- 
blet become possible, Fig. [TJc). We note, however, that 
these features are clearly distinguishable only if E T. 

Current induced magnetic switching. — In order to con- 
trol the adatom/SMM's spin state by a spin-polarized 
current, at least one of electrodes has to be mag- 
netic |12[fT3"]. Here, we choose it to be the tip, Pt ^ and 
P s = 0. Consider first the case of E = 0. When \V\ < T. 
the system's spin can still fluctuate indirectly between 
the states \ ±S) as a result of stepwise 'overbarrier' transi- 
tions. However, due to imbalance of the spin-up and spin- 
down electron tunneling processes, the systems 's spin be- 
comes locked in one of the ground states |±5) [depending 
on the bias polarization] for \V\ > T, Fig[2](d). In con- 
sequence, the conductance is determined in the GMR- 
fashion by the relative orientation of the tip's magnetiza- 
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tion and the adatom's/SMM's spin, dl/dV cx 2aP t (S z ), 
Fig [2jc) . This behavior stems from interference of direct 
tunneling and spin conserving part of tunneling associ- 
ated with exchange interaction between tunneling elec- 
trons and magnetic core of an adatom/SMM, and also 
leads to the asymmetry in conductance with respect to 
bias reversal, see Fig. [2^c). Inelastic transport processes 
activated at |V| ps Vthr enhance conductance [see the 
curve for E/D = in Fig[2jc)], but slightly reduce the 
absolute value of the average spin S z component [see the 
curve for E/D = in Fig^d)]. 

The situation changes qualitatively for E ^ 0. For 
low voltages, the spin oscillates between \x±s), like 
in the nonmagnetic case. However, due to the spin- 
dependence of tunneling processes 7 x _ sXs ^ 7 X sx-s' an d 
the spin generally resides longer in one of the two ground 
states At some voltage, transition rate to the ex- 
cited state surpasses the transition rate between the two 
ground states and the spin becomes locked in one of the 
two states \x±s)- This results in an additional peak in 
dl/dV, Fig. [2^c). Since the transition rate 7 X _ SXS in- 
creases with E, position of this peak moves towards larger 
voltage with augmenting E. Further rise in voltage leads 
to saturation of the conductance, and the saturated value 
only weakly depends on E. It's worth of note, however, 
that the strong mixing of spin states for large transverse 
anisotropy prevents the system's spin from aligning along 
the easy axis. Surprisingly, in the case of spin-polarized 
transport, the interplay of the underbarrier relaxation 
process introduced by the transverse magnetic anisotropy 
imposes the voltage barrier for switching the system's 
spin. Such a behavior doesn't take place in systems ex- 
hibiting purely uniaxial magnetic anisotropy. 

When an external magnetic field is applied, one can 
achieve degeneration of different states belonging to ei- 
ther of the two decoupled manifolds, see the right side 
of Fig. [TJc). This degeneration, though, is not complete 
due to level repulsion, as shown in Fig. |3je). Anyway, at 
the fields corresponding to the dashed lines in Fig. [3|e), 
quantum tunneling of magnetization takes place, which 
results in transitions between the two degenerate states. 
These transitions are clearly seen as resonant peaks in 
the average values of the molecules's spin, Fig. [3ja,c,d), 
and current, Fig.[3jb). The resonant character of QTM is 
even more evident in the field dependence of dI/dB z , see 
Fig. |3jc). Interestingly, the resonances due to QTM can 
be observed only when (at least) one of the electrodes is 
ferromagnetic, as follows from Fig.^d). 

Conclusions. — In this Letter we have considered the 
influence of transverse magnetic anisotropy on spin polar- 
ized transport through magnetic adatoms/molecules, and 
in particular on the current-induced spin switching. First, 
we have demonstrated that mixing of states by trans- 
verse anisotropy leads to a decrease in conductance in the 
elastic transport regime (low voltage regime) and to ap- 
pearing of the peaks at voltages where the system's spin 
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Figure 3. (Color online) (a) The average value of the spin's 
zth component and (b) the charge current shown as functions 
of an external magnetic field applied along the system's easy 
axis, B z , and bias voltage for D = 50 /xeV, E/D — 0.1 and 
Pt = 0.5. (c) Dark bold line represents a cross-section of (a) 
for indicated bias voltage [dashed lines correspond here to the 
case of E — 0], while bright bold line corresponds to d//d£> z , 
shown here in arbitrary units, (d) Average spin for V = l mV, 
E/D = 0.1, and for several values of the tip polarization Pt. 
(e) Dependence of the system's energy spectrum on magnetic 
field B z . Other parameters as in Fig. [I] 



becomes locked in one of the two ground states. When 
an external magnetic field is applied, the phenomenon 
of quantum tunneling of magnetization, which occurs at 
some resonant values of the magnetic field, results in res- 
onant peaks in the average value of the molecules's spin 
and in the charge current. It is worth emphasizing, how- 
ever, that these effects can be observed only when the tip 
(and/or substrate) is ferromagnetic. Thus, spin-polarized 
transport spectroscopy may prove a useful experimental 
tool for studying the phenomenon of QTM. 

It has been also shown that the conductance gen- 
erally depends on the relative orientation of the av- 
erage adatom's/SMM's spin and electrode's magnetic 
moment. This dependence stems from the interference 
of direct tunneling and spin conserving tunneling con- 
nected with exchange interaction between the electrons 
and adatom/SMM spin, and resembles the giant mag- 
netoresistance effect in magnetic multilayers. It is also 
responsible for a significant asymmetry of transport char- 
acteristics with respect to the bias reversal and can be 
used to control spin switching of the adatom's/SMM's 
spin. 
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